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ABSTRACT 



The use of boiling heat transfer in the liquid immersion 
cooling of electronic components has always been hampered by 
the excessive superheat necessary to initiate nucleation in 
the fluorinated hydrocarbons used as dielectric cooling 
fluids. In an attempt to overcome some of these difficulties, 
an experimental study of the effects of nucleate pooling 
boiling, on the heat transfer from surface near the boiling 
surface was conducted. An experimental chamber was constructed 
which had a column of four horizontal wires spaced 2.5 cm 
vertically from each other. The lowest wire was progressively 
heated from the natural convection region through the nucleate 
boiling region. A study was made of the effects of the boiling 
wake from the lowest wire on heat transfer from the upper 
wires. Under certain conditions heat transfer enhancements of 
up to 30% were obtained. 
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I . INTRODUCTION 



A. BACKGROUND 

Research on pool boiling in electronic cooling systems has 
focused on three primary areas: (1) reducing the temperature 
excursion at incipient boiling, (2) reducing wall superheat 
during nucleate boiling and (3) enhancing critical heat flux 
(Mudawar, 1990) . Since the heat transfer potential of nucleate 
boiling is well known, this regime was one natural selection 
for examination as a possible solution to this cooling 
problem. Nishikawa (1980) notes that there are two primary 
methods used to promote nucleate boiling. The first is to 
treat the surface in a manner that reduces its wettability. 
The second, more promising one, is to manufacture a surface 
with cavities which trap and hold vapor keeping the nucleation 
site active. 

In recent years the drive toward higher heat flux and 
power density in electronic equipment has turned attention to 
the use of dielectric liquids. The 3M Corporation has 
commercially produced its "Flourinert" series of high 
dielectric, inert electronic cooling liquids. In contrast to 
water and other conventional coolants these liquids possess 
relatively low thermal conductivities and boiling points and 
small heats of vaporization. They demonstrate high wettability 
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on most surfaces. Their surface tensions are the lowest of any 
known liquids, resulting in near-zero wetting angles. This is 
significant because it eliminates many potential nucleation 
sites. High dielectric strength and high resistivity are some 
of the characteristics that make these liquids well suited for 
electronic applications. Especially, FC-72 is known as 
chemically stable and leaves no residue during boiling. 

However, when using dielectric fluids, incipience 
superheat excursion is pronounced which leads to very high 
wall temperatures. A firm understanding of boiling incipience 
is required in order to implement effective cooling. Boiling 
can produce very large-scale heat transfer coefficients but it 
also introduces some problems: (i) Boiling restricts the 
physical design of the system, (ii) A high degree of superheat 
may be required if the surface is very smooth in order for 
nucleation to occur. Due to the high local heat flux the 
occurrence of periodic boiling may result in local turbulence, 
even damage of electronic components. (iii) Impurities 
concentration may result at the nucleation site. 

Subcooling offers many advantages over saturated boiling 
for cooling large arrays of computer chips due to the 
decreased diameter and rapid collapse of departing bubbles. 

In this study heat transfer enhancement due to bubble 
pumping in FC-72 was investigated. It is hoped that they will 
serve to identify some trends and provide a basis for further 
research on actual components and integrated circuits. 
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B . PREVIOUS WORK 



Hsu (1962) proposed a method for determining the range of 
active nucleation sites as a function of temperature or heat 
flux. This model revealed that the maximum and minimum sizes 
of active cavities are functions of subcooling, system 
pressure, physical properties and the thickness of the 
superheated liquid layer. 

Bergles and Chu (1979) compared the nucleate pool boiling 
heat transfer characteristics of three copper Union Carbide 
High Flux test sections to a plain copper tube in distilled 
water and R-113. The experiments were conducted with three 
different treatments to the test surface prior to collecting 
data: subcooling; aging by preboiling in the pool; and heating 
the surface in air to remove all liquid in the pores. This 
research indicated: 

1. Heat transfer coefficients for the porous surfaces 
were four to ten times higher than for the plain 
surfaces . 

2. There was a significant temperature overshoot 
prior to the initiation of boiling and it was not 
very sensitive to aging or the power increment 
changes in water. These overshoots were sensitive 
to aging, initial subcooling and power increment 
changes in R-113. 

Bergles et al. (1968), conducted nucleate pool boiling 
studies in water, R-113 and FC-78, using a 304 stainless steel 
tube. As heat fluxes were increased, the heat transfer 
coefficient followed the convective curve until a high 
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superheat was attained, and then a distinct increase in the 
heat transfer coefficient was noted as nucleation occurred. 
This was followed by an abrupt drop in all temperatures by as 
much as 27°F. It was also noted that the inception of 
nucleation could be triggered by vibrating the test surface. 
The conclusions drawn were: 

The phenomenon of temperature overshoot hysteresis in 
ordinary liquids is due to two causes: 

1. The existence of metastable bubbles which are 
triggered only at sufficiently high disturbance 
levels , 

2. The deactivation of larger cavities by displacing 
the vapor by liquid during subcooling. 

Many of the experimental reports published prior to 1986 
are reviewed by Bar-Cohen and Simon (1988) . They noted 
significant variations in the incipience superheat from one 
report to the next. The authors discussed possible mechanisms 
for delayed nucleation. Recently non-repeatable and unsteady 
behavior of boiling incipience for a highly-wetting dielectric 
liquid (R-113) was documented by You et al. (1989) . 

The bubble behavior from the nucleate boiling regime to 
the film boiling regime has not been clarified yet, except for 
liquid helium. The transient boiling experiment of liquid 
nitrogen at stepwise heat generation in a platinum wire was 
carried out by Kunito Okuyama (1989) . 
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McAdams et al. (1949), reported the effect of dissolved 
air on forced-convection subcooled nucleate boiling of water 
where the dissolved air reduced the heat flux at boiling 
initiation . 

Torilcai et al. (1970), used platinum wires immersed in 
water to measure dissolved gas effect on the incipience wall 
superheat under reduced pressure conditions (0.05 - 1.0 Bar) . 

Surface enhancement is a very popular tool for improving 
pool boiling performance. Marto and Lepere (1982) found three 
commercially available enhanced boiling surfaces which shifted 
the boiling curve of FC-72 toward lower superheats compared to 
a plain copper surface. In a study by Berenson (1962), 
nucleate boiling heat transfer coefficients varied by 600% due 
to variation in surface finish. He emphasized the importance 
of the surface roughness. Recent investigations on surface 
effects in pool boiling by Chowdhury and Winterton (1985) show 
that heat transfer appears to be unaffected by contact angle. 

The effects of flow velocity and subcooling on critical 
heat flux from a simulated electronic chip attached to the 
wall of a vertical rectangular channel was determined by 
Mudawar and Madox (1988) . They concluded that heat flux was 
significantly enhanced by increasing the degree of subcooling. 

The solubility of gas in highly-wetting dielectric fluids 
is of special concern. These fluids can contain up to 25 times 
more dissolved air by volume than can water: 48% air by volume 
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can be dissolved in FC-72 compared to 1.9% air by volume in 
water (3M Product Manual, 1987) , 

Experimental results presented by You, Simon and Cohen 
(1990) indicate that dissolved air plays an important role at 
incipience only when a very high dissolved air content is 
maintained. The surface tension value for FC-72 was estimated 
to have a 2% reduction due to high dissolved air content. 

S.M. You et al. (1989) conducted an experiment in 
saturated R-113 fluid under atmospheric pressure with two 
types of wires. One is 0.13 mm diameter chromel wire (90% 
nickel) and the second is a thin platinum film on a quartz 
rod. They examined many pits of various sizes using a scanning 
electron microscope (SEM) . 

Bar-Cohen (1990) conducted experiments providing pool 
boiling characteristics for a number of commercially available 
inert liquids (FC-43, FC-72, FC-75, FC-77, FC-84, FC-87) . Data 
were obtained with a classical 0.25 mm platinum wire heating 
element. He concluded that: 

1. Nucleate boiling characteristics of these fluids 
depend on the heating surface, 

2. Saturated pool boiling curves were almost 
identical with one another on a particular heater, 

3. Superheat excursions at boiling incipience appear 
to always occur, but their extent can vary widely. 
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C . OBJECTIVES 



The objective of this work was to investigate the effects 
of the boiling wake from a heated wire on the heat transfer 
from wires placed within the wake. This work was carried out 
using the dielectric liquid FC-72. 
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II. EXPERIMENTAL DESIGN 



A. DESCRIPTION OF COMPONENTS 

Figure (1) shows the system configuration and identifies 
all major components. Figure (2) is a photograph of the 
apparatus . 

The chamber consisted of a plexiglass box with an aluminum 
coverplate fitted with a rubber 0-ring. This vessel has side 
wall heaters to allow the liquid to be preboiled for degassing 
before the run, and to maintain the liquid at saturation 
temperature throughout the run. The vapor was condensed and 
returned to the chamber by gravity from the aluminum cover 
plate condenser equipped with thermoelectric coolers. The 
following is a description of the test apparatus. 

1 . Test Chamber 

The test chamber was constructed of plexi-glass of 1/2 
inch thickness and had inside dimensions of 6 inches height, 
6 inches length and 2 inches width. An isometric view of the 
chamber is presented in Figure (3) . Side walls are glued to 
each other. 

It has a 1 inch diameter hole and 1 inch protruding 
pipe and a copper elbow tube in the back face which is the 
pathway for the power and thermocouple wires. The chamber also 
allows the test section to slide in and out for the 
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experiment. Cross-sectional and schematic drawings of the test 
section in the chamber are presented in Figures (4) and (5) 
respectively. The chamber also contained two small (1/4 in) 
pressure equalization outlets, located at the top-right and 
the bottom-left of the front vertical surface. It is a simply 
a closed tank with an aluminum cover plate. 

The test chamber is attached to a wood base in order 
to give better arrangement for cables and connections. 

2 . Thermocouples 

Thermocouples, 0.01 in (0.254 mm) in diameter, were 
used to take measurements. The thermocouples for the plate 
condensers (thermoelectric coolers which are located on top of 
the aluminum plate) were bonded in grooves by placing a small 
drop of Eastman 910 adhesive on the thermocouple bead. The two 
thermocouples then were placed in the designated line groove. 
Using a weight, pressure was applied until the adhesive set. 
After allowing for three hours cure time, the remainder of the 
groove was filled with Omega Bond 101 Epoxy and smoothed to 
the level of the aluminum plate inside face. The Omega Bond 
101 Epoxy experienced little shrink after curing. 

3 . Heaters 

The side wall heaters were glued on either side walls 
of the chamber. The actual size of the heaters were 125 x 11 
mm and they were subjected to 5 V, 2.5 amperes. The power load 
depended on the combination of the wires being heated. 
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4 . Heat Exchangers 

The three thermoelectric heat exchangers measure 1.56 
X 1.56 inches and were placed on top of the aluminum cover 
plate. The power supply for these were determined by the 
related curves. The property curves and technical data are 
presented in Figures (6) and (7) respectively. 

5 . Al\imin\am Cover Plate 

An aluminum plate of 1/2 inch thickness was used as a 
cover plate for the test chamber. The inside face of the plate 
had two thermocouple beads to monitor the condenser 
temperature. Three thermoelectric condensers were placed on 
top of the aluminum plate. A thin groove was cut around the 
edges of the inside face for a 0-ring gasket. The plate was 
screwed over the test chamber with eight stainless steel 
screws of 1/16 inch diameter and compressed the 0-ring to form 
a watertight seal. Figure (8) shows the cover from the top 
with the condensers on it . 

B. INSTRUMENTATION 

Separate voltage taps across the heating element and 
across the precision resistor are used to compute heater 
resistance and heating element power. The heating element 
temperature is computed from its resistance via temperature 
calibration curve obtained prior to the experiment. 
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1 . Power to the Heaters 



Each heater was run in series with a precision 
resistor that was measured to have a 0.1 ohm resistance. All 
four heaters were in parallel with a Hewlett Packard model 
number 6289A, 0-40 volts, 0-1.5 amperes, DC power supply. Both 
the voltage across the precision resistor and the chromel wire 
were measured independently. The current to each heater was 
calculated by the product of the heater voltage and heater 
current. Both voltages were measured by a Hewlett Packard 
model 3852S data acquisition system containing a Hewlett 
Packard model 44701A digital voltmeter, all controlled by a 
Hewlett Packard model 300 computer. 

2 . Temperature Measurement 

The thermocouples for. the bulk fluid, ambient, ice 
reference point and coolers inside plate temperatures were 
referenced individually to an electronic Ice Point Reference 
as seen in Figure (9) . Each reference thermocouple was 
connected such that its constantan lead was connected to a 
constantan lead of the measurement thermocouple. The copper 
leads of each measurement and reference thermocouple were 
connected directly to a Hewlett Packard model 44705A 20 
channel relay multiplexer and inserted into the data 
acquisition system. The data acquisition system then measures 
the ice referenced voltage of each thermocouple. The voltages 
were converted directly into temperatures in the controlling 
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computer program by using a fourth order polynomial, fit to 
the thermocouple manufacturer's calibration data for copper- 
constantan thermocouples (Hazard, 1986) . The locations of the 
thermocouples in the chamber are presented in Figure (10) . 
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Figure 1. System Configuration. 




Figure 2 . Photograph of Test Apparatus . 
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Figure 3. Isometric View of the Test Chamber. 
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Figur© 4. Cross Sectional View From th© Top of th© Test Section. 




Key: 

1. Washers 

2. Chromel Wire 

3. Screw 

4 . Nut 

5. Power Line 



Figure 5. Schematic Drawing of Test Section in Vessel. 
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Figure 7. Technical Data for Thermoelectric Coolers. 
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Figure 8. Cover Plate with Thermoelectric Coolers. 
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Figure 9 . Thermocouple Connection Schematic . 




Figure 10. The Locations of the Thermocouples in the Chamber. 
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III. EXPERIMENTAL PROCEDURE 



A. NORMAL OPERATION 

In carrying out the experiments, the procedure was as 
follows ; 

1. The ice point reference device was turned on. A mercury 
thermometer inserted into the ice base ensured that the ice 
bath is at 0 ± 0.05°C within one hour. 

2. The data acquisition system was then turned on. The 
internal voltmeter in the data acquisition system requires one 
hour warm-up period. 

3. The computer and printer were turned on. The data 
acquisition program contained -in Appendix D was loaded. 

4. The HP 6289A, 0-40 volt, DC voltage supply was set to 
zero and energized. The current limiter on the 40 volt power 
supply was set to the maximum current level. 

A sampling of all thermocouple temperatures was performed 
and printed. The power was increased incrementally and the 
output from the power calculation program was checked. 

The bulk fluid temperature as given by the wire resistance 
must match the value given by the bulk fluid thermocouples. 
When the results showed that the temperature stratification to 
be less than 0.10°C the experiment was ready to begin. The 
power to the test section is then increased by controlling the 
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test section heater voltage. The voltage was increased in 0.2 
volt increments until the initiation of boiling occurred, at 
which time 0.1 volt increments were made until the maximum 
attainable heat flux was reached. This occurred at 5.7 volt 
for the power supply, making use of the HP 6289A, 0-40 volt, 
0-1.5 Amperes DC voltage supply up to 8.1 volts, which 
corresponds to the maximum value. The power supply has a reset 
function to protect the device itself. The power is then 
decreased in the same fashion. At each power setting the 
system was permitted to stabilize for a few minutes and the 
following data were recorded: The heater voltage, precision 
resistor voltage, two thermocouple indications of 
thermoelectric cooler temperatures, and three thermocouple 
readings of pool temperatures. This procedure was repeated for 
every set up. 

B . DATA REDUCTION 

The test wire temperature was found using the calibration 
formula between thermal resistance and temperature value of 
the chromel wire. 

R = R^(1+eAT) 

where. 

At = T^ - Tj, 

is the wall temperature 
Tfc is the bulk liquid temperature 
e is the slope of the calibration curve. 
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Details of the calibration slope value are presented in 
Appendix C. This study utilizes the value of 0.0004 C/ohm 
offered by Kaye (1976) . 

The starting value of the wire resistance is accepted as 
the in the previous formula. 

The heat flux is calculated from the following formula; 

Flux = Power /A 

where, A = k D L is the total wetted surface area, and D and 
L are diameter and length of the wire respectively. 

The bulk temperature of the liquid was measured by taking 
the arithmetic average of the three thermocouple readings. 

Ttui, = ( T, + + T3 ) / 3 

The Rayleigh number based on diameter is given as: 

^ gp ( T,-r,) 

° V a 

where, 

P = 

^ TpTcpT 

All the properties in the above equations are calculated 
at the bulk temperature. 

Determination of Nusselt number is: 




where h is the heat transfer coefficient and was determined 
from 
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^ “^bulk ) 

where, 

= wire surface temperature (°C) 
h = heat transfer coefficient (W/m^ • C) 
g" = heat flux (W/m^) 

Churchill and Chu (1975) have recommended a single 
correlation formula for the natural convection from long 
horizontal cylinders. 



Nu^={ 0.60 + 



0.387 Ray^ 



[ 1 + ( 0.559/Pr ) 



T7TTTT777 



It is valid for a wide Rayleigh number range 

10'^ < Rad < 10^^ 

where the Prandtl number (Pr) is defined as 
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IV. RESULTS AND DISCUSSION 



Data were taken and the related graphs can be grouped for 
discussion purposes as: 

1. individually powered wires 

2. pumping effect runs 

A. INDIVIDUALLY POWERED 

Before investigating the pumping effect of the bubbles, 
individual boiling experiments on the wire heaters were 
conducted. 

Figure (11) shows wire 1 (the lowest wire) where distinct 
single phase and nucleate boiling regimes are observed. When 
the heat flux was increased in small steps, the single phase 
curve reaches 17°C superheat. At about 20°C an abrupt 
temperature decrease occurs as nucleate boiling begins. 

Figures (11) and (12) are essentially the same. The 
coordinates are heat flux versus surface temperature of the 
chromel wire. Figure (12) is on the log-log scale while Figure 
(11) is semilog graph. 

Figure (13) is the single run boiling curve of heater 
number 3. It was the only one to show overshoot of less that 
7°C. 
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Figure (14) for wire number 1 shows the same type of 
behavior as Figure (11) . In this case data were taken at 
different times in order to demonstrate the repeatability of 
the experiments. Data points for two different runs did not 
coincide, but that is because the bulk fluid temperature is 
not the same for all experiments, but depends on the ambient 
temperature of the room. 

Figure (15) indicates all individually powered heaters on 
the same graph. The difference in the overshoot may be because 
of the surface geometry of heater number 3. 

The large aspect ratio of the heaters resulted in minimal 
conduction losses to the supports and yielded a uniform 
temperature along the length. 

The boiling curves for all four wires are practically 
superimposable, even at the departure from nucleate boiling. 
They differ only in the heat flux at which boiling incipience 
occurred and at which full nucleate boiling began. 

Figure (16) shows the best fit to the natural convection 
portion of the experimental data for wire number 1, and is 
compared with Churchill (1975) correlation values. There is a 
constant offset from the correlation values because of the 
wall efffect except for the low Rayleigh numbers of that run 
which also shows considerable scatter in the data. 

The saturation temperature of the FC-72 is about 56°C. 
^sac “ “^bulk is the subcooling. Figure (17) shows the effects of 
subcooling. At boiling incipience, a liquid mass near the 
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heated surface is slightly superheated, even though the bulk 
fluid is subcooled. Each run began with zero heating. The heat 
flux value was increased incrementally, continuing up the 
single-phase natural convection line until boiling was 
initiated. The figure indicates that as the subcooling 
increases the curve shifts to the left, displaying significant 
increases in the overshoot and heat fluxes. With bulk liquid 
at 15°C, it shows large amount of superheat, probably a bubble 
begins to form on the wire and as soon as the bubble grows a 
bit it encounters subcooled liquid and collapses. The liquid 
is so cold that it doesn't allow the bubble to grow. As the 
bulk temperature of the liquid gets warmer, it allows the 
bubbles to grow. Subcooling results in smaller bubbles, which 
are less likely to spread to neighboring nucleation sites than 
the larger bubbles associated with saturated boiling. 

B. PUMPING EFFECT 

For this type of experiment, heater number 1 was used as 
the pumper. 

Wire number 4 was fixed at various heat flux values, while 
heater number 1 was increased in a stepwise manner. 

Figures (18) and (19) show the pumping effect of the 
bubbles generated by wire number 1 on the temperature of wire 
number 4. Clearly Figure (19) gives better understanding by 
way of marked corresponding points. As the temperature comes 
closer to the temperature at which the boiling incipience 
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starts, wire number 4 is little effected. Just after wire 
number 1 is subjected to overshoot it also has similar a 
experience, i.e., a sharp jump . Actually, it was the general 
pattern for most of the other runs. This suggests that boiling 
can be started by the pumping effects of the rising bubbles 
which trigger the heater located above it. The characteristics 
of the bubbles leaving the lower heater are uncertain. They 
may be the wake of the buoyancy-induced motion, especially 
when the upper wire was fixed at a high heat value. A similar 
trend can be observed from Figure (20) . One thing to notice is 
that heater number 4 is fixed at a heat flux that corresponds 
to a starting value of surface temperature of SOX for this 
run. This is below the temperature that meets the overshoot. 
Less jump in the temperature is associated with the lower 
fixed heat flux value. If we compare it with Figure (19), 
which was at higher heat value than this one. Figure (21) 
provides a detailed look at Figure (20) . 

Figure (22) shows wire number 3 fixed at 37,490 W/m^. It 
has the points marked in chronological order. The 3rd wire 
heater is located closer to the first wire than is the fourth 
wire, so it allows most of the bubbles to reach the upper 
heater with collapsing due to surrounding subcooled liquid FC- 
72. 

Figure (23) shows another combination of heater number 3 
with heater number 1. Heater number 3 is kept at 20,450 W/m^. 
This was very low heat flux value that the triggering effect 
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did not show itself. Bubbles still caused some decrease in the 



temperature of the upper wire. 

Figures (24a) and (24b) show very similar results. In this 
case wire number 2 is Jcept at 20,530 W/m^. Since wire number 
3 is the closest to wire number 1. Being so close helped it to 
lose some heat corresponding to a decrease in its temperature. 

Figures (25a) and (25b) show no effect from the bubbles 
from below. Here heater number 2 is kept at a very high heat 
flux of 212,700 W/m^, which is already boiling. It is not 
sensitive any more to the bubbles produced by heater number 1. 
Consequently no jumps or apparent decrease on the surface 
temperature of the number 4 wire are seen. 

Figure (26) shows similar trend as the early graphs. In 
this case the bubbles are produced by wire number 2 and the 
effects on wire number 3 are presented. Heater number 3 had 
the temperature decrease as the pumper had overshoot 
experience . 



31 



HEAT FLUX (W/M 2) 



~o 

E 



































■ 










■ 

■ 

■ 


















1 1 




1 






i ■ — t 

1 


di 








^ ■ 












■ 










■ 












m 

■ 










■ 

■ 










■ I 




1 




) 1 1 








1 


1 


1 








■ 1 


1 


1 




















1 












■ 


















1 




LEGEND 

■ = HEATER WIRE #1 - 
















1 










1 










1 















23.0 34.0 45.0 56.0 67.0 78.0 89.0 



TSURFACE (C) 

Figure 11, Wire Number 1 Individually Powered, Semi-log. 
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Figure 12. Number 1 Heat Flux vs. Surface Temperature, 

Log-log . 
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Figure 13. Number 3 Shows Less Overshoot. 
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Figure 14. Number 1 at Different Times. 
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Figure 16. Best Fit Experimental Data Compared with 
Churchill (1975) Correlation. 
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Figure 17. Subcooling. 
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WIRE #4 AT CONS.HEAT FLUX 
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Figure 18. Pumping Effect with Wire Number 4 at 42,800 W/m^. 
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Figure 19. Pumping Effect with Wire Nximber 4 
(Detailed Figure) . 
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WIRE #4 AT CONST. HEAT FLUX 
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Figure 20. Wire Number 4 at Constant Heat Flux (38, 600 W/m^) . 
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Figure 21. Detail of Figure 20. 
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Figure 22 . Wire Ntunber 3 at Constant Heat Flux Value of 

37,490 W/m^ 



43 



HEAT FLUX (W/M 2) 

■4 



WIRE #3 AT CONS. HEAT FLUX 

































■ ■ 










■ 


■ 

■ 










j 


f “ 










■ 


























■ ■ 


r 










■ 










■ ■ 










■ 

'mm 


■ 








■ 
















1 

1 


t 1 

1 i 










1 


1 


1 








1 


1 


1 








1 


1 






— ■ 






1 


















■ 


















LEGEND 

■ = #1 INCREASING 
o = WIRE #3 FIXED 





































































^ 23.0 34.0 45.0 56.0 67.0 78.0 89.0 

TSURFACE (C) 




Figure 23. Wire Number 3 at 20,430 W/m^, 
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Figure 24a. Wire Number 2 at 20,530 W/m^. 
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Figure 24b. Detail of Figure 24a. 
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Figure 25a. Wire Number 2 Kept at 212,700 W/m^. 
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WIRE #2 AT CONS.HEAT FLUX 




Figure 25b. Detail of Figure 25a. 
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Figuxe 26. Wire Number 2 with Wire Number 3 at 
Constant Heat Flux. 
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V. CONCLUSIONS 



The figures give the idea of the pumping effects. When the 
pumping heater is powered progressively, the upper heater will 
react as indicated in the following. 

• If the upper wire is located close to the lower 
pumper, it will meet most of the bubbles. It will lose some 
heat depending on how close it is. 

• The initial heat flux value of the upper heater is 
also important. The amount of temperature jump which 
simultaneously occurs when the pumper has overshoot depends on 
how warm the upper wire is. 

• If the upper heater already has bubbles, this means it 
had incipience before, and the bubbles will not help it lose 
more energy. 
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VI. RECOMMENDATIONS 



1. Use platinum wires rather than chromel wires to get 
better temperature measurement . 

2. Before the start of each run, boil the pool for a 
sufficient time to degas the liquid. 

3. Modify the software for taking data by adding a WAIT 
command to obtain the time function. 

4. An addition of a foil type heater placed horizontally 
at the bottom of the chamber may help create a uniform 
temperature throughout the bulk fluid. 

5. An extension of tygon tubing of about one foot length 
would take care of overflow and hydrostatic head determined by 
the liquid level inside the filling tank. 

6. The present study deals with only subcooling effects. 
Saturated pool boiling should also be studied. In order to 
have constant temperature in the bulk fluid, the thermo- 
electric condensers and the bulk heaters can be run 
simultaneously . 
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APPENDIX A. SAMPLE CALCULATIONS 



A. DETERMINATION OF INPUT POWER 

For the sample calculations, the output values of wire 
number 1 were used: 



Power ( 1 ) = ^ ( 11 ) 

where volt (12) and volt (11) are the voltages across the 
precision resistor and wire respectively. 



0.0048 

0.1 



1 . 988 ) = 0.96 Watts 



B. CONVERSION FROM VOLTAGES TO TEMPERATURES (®C) FOR 
THERMOCOUPLES 



The HP 3852 A Data Acquisition System channels 1-9 record 
the voltage readings in millivolts. The following polynomial 
formula is for copper constantan thermocouples. 



T(°C) = 0.100860 91 + ( 25727.9 * emf) - ( 7 6734 3.5 * emf^ ) 
+ ( 7802.5596 * emf^ ) - ( 9247486589.6 * + 

( 6. 98E11 * emf^ ) - ( 2.66E13 * emf^ ) + ( 3.SAE14 * emf^ ) 

Using thermocouple 1 

emf = 0.000779 volt T=19.72°C 



C. CALCULATION OF WIRE TEMPERATURE 

The calibration curve for the chromel wire (Kaye, 1976) 
yields the slope e of 4.10”^ ohm/°C. 
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i^ = «„[i+er-Tj 

Here and R are the resistances of the wire at the reference 
point and at any point to be measured R^ and ohm(l) are the 
corresponding variables used in the software. 

T,(l) = ( ( (ohm(l) /RJ - 1) / Kaye ) + 

where, 

T^(l) = wire temperature 

Kaye = 0.0004 ohm/°C constant slope 

To = wire reference temperature (°C) . 

T^(l) = (( (4.1009/4.075) - 1) / 4E-4) + 19.64 

= 35.5°C 

D. CALCULATION OF THE AVERAGE BULK TEMPERATURE 

The average temperature of the liquid was simply the 
arithmetic average of the three thermocouple values immersed 
in the bulk fluid. 

= (T(l) + T (2) + T (3)) / 3 
= (19.9 + 19.54 + 19.6) / 3 = 19.68°C 

E. DETERMINATION OF HEAT FLUX 

HFLUX ( 1 ) 

AREA 

where, 

AREA = K * D * L ^ (3.1415926) D L 
D, L = Diameter and length of the wire. 

AREA = (3.1415926) (0.005) (0.0254^) (3) = 3.04E-5 

HFLUX = = 774 60 W/ 

3.U4E-5 
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F. FLUID PROPERTIES DETERMINATION 



1. Thermal Conductivity k(W/m°C) 

From Figure 5 of the 3M Corporation Flourinert Product 
Manual, the thermal conductivity coefficient curves have been 
determined to be 

fC= ( 0.6033 - 0.00115 ♦ )/ 10 j = 0 ♦ 1 — 

* Cm°C mK 

at Ti, = 50°C it yields 0.0545 W/m. 

2 . Liquid Density p (kg/m^) 

Using the expression in Table 4B and constants 
presented in Table 4C of the Product Manual yields 

p = (1.740 - 0.00261 T^) * 1000 

Tb film temperature must be in units of °C. At T^ = 50°C this 
yields p= 1609.5 kg/m^. 

3. Kinematic Viscosity v (m^/s) 

From Figure (3) and determining a 4th order curve fit 

yields 

V = (0.54 - 0.0116 * Tt + 0.0002085 * tJ) * lE-6 
at Ti, = 50°C this yields v = 4.8E-7 mVs. 

4 . Density of Vapor p, (kg/m^) 

From Table 13A of the Product Manual, values for vapor 
density can be taken. 

For the range of 20-80°C the data can be approximated 
as a fourth order polynomial . 

= (246 - 10.93 * Ti, + 0.213 * t/ - 
0.00193 * + 6.55E-6 * T/J / 338 
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= 50°C it yields 

Pv = 10.67 (kg/m^ 



where 



1 



Pv ‘ 



5. Specific Heat C^, (J/kg°C) 

For all Flourinert Electrochemicals (Figure 4) 

Cp = (0.241111 + 3.70374E-4 * TJ * 4186 
at r = 50°C. Calculation yields 1086.7 J/kg°C. 

6. Thermal Diffusivity a (m^/s) 

a = k/pCp 

at T = 50°C, calculation yields 3.115E-8 m^/s. 

7. Thermal Expansion Coefficient P (1/°C) 

Table 4B in the Product Manual shows 

P = 0.00261 / (1.740 - 0.00261) * TJ 

at T = 50°C. Calculation results in 0.00162 (1/°C. 

8. Prandtl Number (Pr) 

Pr = v/a 

at T=50°C.This yields Pr = 15.4 

9. Calculation of Latent Heat hfg (J/kg) 

Table 13A of the Product Manual shows the values can 
be approximated by a linear equation. 

hfg = (8057 - 22.73 * TJ * 10.79 
where 10.79 is the conversion factor from Cal/Mole to (J/kg). 
At T = 50°C calculation yields 74672 (J/kg) . 
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APPENDIX B. UNCERTAINTY ANALYSIS 



A. UNCERTAINTY IN AREA 

cOfl = .0001 inches 

0)^ = ±0.010 inches 

A = n D L 



where : 
D 
L 
A 



dA 



Wire diameter 
Length of heater wire 
(3.141593 ♦ 0.005 ♦ 3) 
0.0471 in^ 

3.04 * 10'^ 

or L) 



dA 

~Jl 



{n D) 



Therefore, 



^ ) ^ + ( nDdi^) ^ ] 



1/2 



0). 



^ V 

ay. 

K~ 7 



V 

j j 



1/2 
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/ 0.0001 


0.020 


[\ 0.0 05 / \ 


~^l J 



^ = 0.021 or 2.1% 

A 

(0 = 6.38 X 10-® in" = 9.89 x lO'^ m" 

B. UNCERTAINTY IN POWER 

= ±0.05 Ampere Taken from (Hazard, 1986) 

calculation) 

= ±0.5 volt 

= VI 

= I 
= 



0)i 

0 

dQ 

w 

30 

■3T 



Wo = [ ( ICOJ" + ( VCOj)" 



CO 



0 



0 



0.05 \" 

774^1 



/ 0.001 

\ 0.018 



11/2 



= 0.059 or 5.9% 
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c. 



UNCERTAINTY IN HEAT FLUX 



<7 



Q 

A 



dg _ 1 

■Jo ~A 



dg ^ _ 2 

■Ja 



cog 




cog 

~g 




Y 

y 





Y 

y 



11/2 



= [ 0.021' + 0.059' 

= 0.062 or 6.2% 
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D. UNCERTAINTY IN TEMPERATURE 



Using thermocouples: 
Variable 

Voltmeter resolution 

Ice reference bath 
Temperature 

Polynomial Temperature 

Conversion 

R 



Uncertainty 
0.02 6°C 

1.0 mv resolution 
0.05°C 

0.00663°C RMS 



1 . 0 % 



WTi, = [ ( ) 2 + ( M + (d^urve 

= [ (0.0025)^ + (0.05)^ + 0.00663 

coTj, = 0.063 or 6.3% 



E. UNCERTAINTY OF WIRE SURFACE TEMPERATURE 



From the calibration formula 

K = [ 1 + e AT ] (1) 

where, 

AT = - Tfc 

£ = Slope of the calibration curve 

Ti, = Bulk fluid temperature. 

From (1) 






/ e 



For: T^ = 44.14°C 
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1 



23.9°C 






Tb = 



ar, 

~Jr 



1 



e = 0.0004 



a _ -R 

■3T^ " 



= 3.88 Q 



ar, 

aa 






R = 3.91 Q. 



cor. 



0^ = ci^o = C 



{ 



^ar, Y 



( CO r. ) 2 + 



Y 

° J 



(C 0 R„) ^ + 



^ar, 






( CO^ ) 2 + 



^ar Y 



^-3T ^ 



( co^ ) ‘ 



cor. 



{ 1 ( 0 . 063 )^ + ( 6 . 49 )^ + 6 . 44 ^ + 2.31 
9.43°C 



F. UNCERTAINTY IN AT 



Where AT = 



dAT 

~wr 



12 % 



.01 



} 1/2 
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dAT 

-drl 



COat 






dAT.. 

-T 

V / 



Y f 
+ 






dAT,^ 

_, — co~ 

^t' ^■ 

V ^ y j 



1/2 





V 

/ 




* 1/2 
Y 

y - 




= 0.21 or 21% 

COat = 0.21 (44.14 - 23.99) 

= 4 . 3°C 



G. UNCERTAINTY IN HEAT TRANSFER COEFFICIENT 



h 



<7 

~Kt 



q is the heat flux 



dh ^ 1 

dg ~Et 



q = 20190 W/m^ 
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dh 



-g 

Ar^ 



AT = 44.14 - 23.99 






g 



1/2 



-■ .^.,Ci).„ 



= 20.15°C 



TT 






( 0 , 

Lv y 



Y 



g 



^AT 

~KT 




( 0.062 ) 



4.3 

Tins 



ni/2 



= 0.222 or 22% 
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APPENDIX C. CALIBRATION OF CHROMEL WIRES 



The chromel wires of about 3 inch length were immersed in 
a Rosemount constant temperature distilled water bath, and 
were compared with a platinum resistance thermometer as a 
standard. The temperature was raised from 30°C to 90°C in 
increments of about five degrees. At each increment, the 
temperature was allowed to stabilize for five minutes, then 
measured using a Newport digital thermistor devise for the 
0.005 inch diameter chromel wire. All recorded data are 
presented in the following table. The corresponding best fit 
line is also presented on the graph. Two sets of data were 
taken at different times. 



Thermistor 

(C) 


T/C 


Resistance 

(Ohm) 


29 


30.2 


4.74 


38.3 


39.7 


4.76 


49.5 


51 


4.78 


61 . 9 


63.4 


4.8 


76.2 


77.2 


4.83 


87.7 


88.4 


4.85 


52.3 


53.8 


4.79 


57.5 


59.1 


4.8 


60.5 
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4.8 


68.3 


69.6 


4.81 


70 . 7 


72 


4.82 


71.3 


72.4 


4.82 


86.2 


86.9 
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Figure Cl. Chromel Wire Calibration. 
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APPENDIX D. SOFTWARE 



Software written in the BASIC language converts voltage 
measurements into temperature values for the bulk fluid, ice 
reference, and coolers inside the thermocouples. It uses the 
a value of 0.0004 Ohm/°C to calculate the surface temperature 
offered by Kaye (1976) . The printed output covers bulk fluid, 
wire surface temperatures and heat flux. 

A FORTRAN program reads the output values created by the 
previous program and calculates liquid properties R^, Nu, and 

^sf • 

The GRAPHER program graphs the data and prints it using a 
HP-7440A Plotter. 

Initial experiment data was fed into EASYPLOT software on 
the Naval Postgraduate School's main frame computer via floppy 
diskettes . 
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n 
12 
1 3 
U 

15 

16 
16 

19 

20 

25 

26 
27 
29 
50 

31 

32 

33 
35 
35 
37 
35 
35 
40 
4 1 

42 

43 

44 

45 

46 

47 

48 
45 

50 

51 

52 
55 

54 

55 
55 
57 

53 
55 
60 
6 1 



OPTION Bn5E ) 

BEEF 

PRINT "INPUT DnTE" 

INPUT TarihS 
PRINT " INPUT TIME " 

INPUT Zan 

IFRINT "SCANNING MEnS" 

REAL Enf ( 7 ) ,T( 7 } ,D( 3 ) ,Uol t ( 1 I *: 1 9 > ,Pouer( 4 ) ,Hfiux(4 ),Ts(4) ,Ohn( 4 ) 
REAL Sumo 

OUTPUT 705i "CONPHEAS OQU , 1 01 - 1 07 ,U5E 0" 

FOR I«1 TO 7 
ENTER 709;£mf(I) 
j PRINT I ; Enf ( I ) 

NEXT I 

OUTPUT 705; "CONFHEA5 DCv , 1 1 i -1 i 9 .USE 0" 

FOR I«1 1 TO 19 
ENTER 705;Uoit(I) 
i PRINT I ;Uoi t C I ) 

NEXT I 

D< 1 )* . 1 009505 i 
D( 2 >*25727 . 9 
0(3 >*-757345 . S 
0(4 >*78025555 
0(5 >*-9247455595 
D( 5 >*6 . 95E-I- 1 I 
0(7 )=-2 , 5BE+ 1 5 
D(5)*3.S4E+14 
FOR 1*1 TO 7 



Suno*0, 

FOR J=1 TO 5 

5umo*5umoTu( J)«Enf(I>*(J-l ) 
NEXT J 



T( I )*3uno 

Tb*( T( 2 >+T( 3 >+T( 1 > )/3 
NEXT I 

PRINT "DATE ;";Tarih$ 



PRINT 

PRINT 



’TIME 
”T"n . 



; ; lan 

T ( I > ; ” C 



UPPER PORTION OF THE FLUID" 



(PRINT USING ”T"j 1 .DD.0?T(1 C UPPER PORTION OF 
PRINT "T";2 ,T(2>;" C LOWER PORTION OF THE FLUID" 

PRINT -T”;3.T(3>r C LOWER PORTION OF THE FLUID" 

PRINT "T";4.T(4>r C TE COOLER INSIDE AL. PLATE" 

PRINT "T";5,T\5);" C TE COOLER INSIDE AL. PLATE" 

PRINT "T";G,T(G>r C AMBIENT (ROOM TEMPERATURE)" 

TEMPERATURE)" 



THE FLUID" 



66 



d2 print "T" ;7 .7(7 ) ; '’C ICE (REFERENCE TEMPERnTURE ) " 

53 ‘NEXT I 

54 PRINT "T5ULK--’;Tb 

55 jpower coicuiatlon 3 »^*»»»^**»*» 

70 Power( 1 )-(UoIt( 12 )/. ) ;*Uoit( m 
80 Power(3 )»< Voii( 17)/. 1 )*\?oit( I'B ) 

51 Ihealer surface area calculations 

52 i♦♦***♦♦♦*1^** 

53 Dia“ .005* . 0254 

54 nl»3*.0254 

35 ni an*( 3 . 1 4 1 5325 ) *0 i a * n 1 
57 !heat flu>: calc. 

98 

99 Hflux( 1 )«Power( 1 )/Aian 

100 Hf lux( 2 )*Fower( 2 )/Aian 
• 01 Hf iux( 3 )=Fouer( 3 ;/Mlan 
iOZ PRINT '"power 1 = ”; FouerC 1 ; 

103 IPRINT *'power2"“ ; Fouer( 2 ) 

104 IPRINT "Dower5“ Power ( j ) 

105 l*-»***^^*->-r*»-*-»^^^***^^^*-»** 

^ 0hn( 1 )*\joi t (11 )/( Uoi t •: 1 2 > / . 1 ) 

1 05 !0hn(2 )=Uoi t ( 1 4 )/( volt( 1 5 )/ . 1 ; 

uhn(3)=Uolt(lE)/i Uoi t ( 1 7 )/ , 1 ; 

111 IPRINT "0HH2* “ ; 0hm( 2 ) 

112 I PRINT "OHM 1 " " ; Chm( 1 ) 

114 IPRINT ■’0HM3=" ;0hn( 3 ) 

1 1 5 Dt=Tb-Ta 

115 ! Ta 15 THE TEMP OF BULK FLUID Tb UIHEN POWER OFF 
117 ! Ro IS THE RE5I5. OF CHROhEL WIRE nCTUALLY IT IS 

115 i THE^Ohn(l) UnLUE WHEN THE SULK FLUID TEMP (Tb-Tanb) AND POWER OFF 
119 I ( POWER OFF MEANS Power ( 1 ;"e . 0054 and it is obtained by setting ) 
1*^® ■ F0wER( 4 ) *0 . 0088 

121 ! R0 1*4. 05 9 R04*3.7931 

Ijj i (UiZ 40 U Supplyist to 00 . 1 0.10 position WHEN 15.4 INSERTED 

124 I Kaye is the value of the calibration slope ReflS.M YOU et all 

125 1 

125 Ta*Tb 

127 Roi=4.0B9 

125 Ro3*3.7851 

123 Kaye*. 0004 

130 R=4 . 05 9* ( 1 *!• , 0004 ^Oi ) 

1 o 1 1 s( 1 ;*( \ ( Ohn\ 1 )/Ro 1 ) - 1 )/Kaye ;*1 ‘Tq 

1j2 is(o;*((( Ohrj( 3 )/ Ro3 )"*! )/Kaye) + Ta 

133 Tf iin*( Tb+Ts( 4 ) )/ 2 

134 PRINT "TSi •" ;Tsi 1 ; . " KPLUX {"jHfluxH) 

135 PRINT "TSS"" ;Ts<3) HPLUX3 ;";Hflux<3> 

13S !*!HINT "T flin ;";Tfiin 

340 End 
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* 



•» 

* 



PARAMETER ( N-ES ) 
INTEGER I 



REAL TS (N) , TB (N) , HF (N) , NU (N) , RA (N) , NUC (N) , TF (N) , TFl (N) , BETA(N) 

,DENF(N),DENU(N),KINV(N),CP(N),K(N) ,HFG(N), ALF(N),PR(N) 

, SURT (N) , RAYD (N) , RAYFB (N) , NUSEX (N) , NUSC (N) , DIA, G, V (N) 

, MU (N) , Y (N) , X (N) , UlS (N) , TSAT, CSF (N) 



READ Tsurface, Heat Flux,Tbulk FROM DATA FILE * 

WRITE Tsurface, Tbulk, Heat Flux, Nusse It , Ray le igh, Nusselt (Chur) * 
Tf i Im , beta, density, viscosity, Cp, conductivity, latent heat, 
Prant 1, surface tension 



OPEN (UNIT=»1,FILE=»’ INPUT. DAT’ , STATUS^’ OLD’ ) 
OPEN (UNIT=2, FILE=»’ OUTPUT. DAT’ , STATUS-’ NEW’ ) 
OPEN (UNIT-3, FILE-’ PROPEl. DAT’ , STATUS-’ NEW’ ) 
OPEN (UNIT-4, FILE-’ PR0PE2. DAT’ , STATUS-’ NEW’ ) 
OPEN (UNIT-5, FILE-’ ROHSE. DAT’ , STATUS-’ NEW’ ) 



DO 15 I = 1,N 

READ (1,*) TS(I),HF(I),TB(I) 

TF(I) = (TS(I)-^TB(I) )/2. 

TFl (I)-TF(I)+273. 15 

BETA ( I ) =t3. 0026 1/(1. 740-0. 0026 1»TB ( I ) ) 
PRINT*, ’BETA IS -’,BETA(I) 

D I A-0. 005*0.0254 

G-s.ai 

TSAT-56. 

WS(I)=TS(I)-TSAT 

DENF ( I ) = ( 1 . 740-0. 0026 1 *T0 ( I ) ) * 1 000. 



U ( I ) = ( 246- 10. 93*TB ( I ) +0. 2 1 3*TB ( I ) **2- 19 . 3E-4*TB ( I ) **3+ 
+ 6. 55E-6*TB( I ) **4) /338. 

DENV (I)-l. /V (I) 

KINV ( I ) - (0. 54-0. 01 16*TB ( I) +0. 0002085*T5 ( I ) **2) *lE-6 
CP ( I ) = ( 0. 24 1 1 1 1 +3. 7037E-4*TB ( I ) ) *4 1 86 . 

K(I)-(-0. 00115*TF(I)+0. 6033) /10. 

HFG ( I ) = ( 8057-22. 73*TB ( I ) ) *10. 79 



ALF ( I ) =K ( I ) / ( DENF ( I ) *CP ( I ) ) 

PR(I)=KINV(I)/ALF(I) 

SURT(I)-(-0. 12*TF(I)+15)/1000. 

MU ( I ) -KINV ( I ) *DENV ( I ) 

C **-*********HM*”****************************************** 

Y(I) = (HF (I) / (MU( I)*HFG( I) ) ) *SQRT(SURT ( I) / (G* ( DENF ( I ) -DENV ( I) ) ) ) 
X ( I ) =CP ( I ) * (TS ( I ) -TB ( I ) ) / (HFG ( I ) *PR ( I ) **1 . 7) 

CSF ( I ) =X ( I ) / ( Y ( I ) **0. 333) 




RAYD ( I ) =G*BETA ( I ) * (TS ( I ) -TB ( I ) ) *DIA**3/ (KINV ( I ) *ALF ( I ) ) 
RAYFB ( I ) =G*BETA ( I ) *HF ( I ) *DIA**4/ (K ( I ) *KINV ( I ) *ALF ( I ) ) 

NUSEX (I)=HF(I)*DIA/ (K ( I ) * (TS ( I ) -TB ( I ) ) ) 

NUSC ( I ) - ( 0. 60+ ( 0. 387*RAYD (I) **(l./6. )/((! + ( 0. 559 /PR (!))♦♦( 
9. /16. ) )**(8. /27. ) ) ) )**2 



WRITE (2,20) 
WRITE (3,30) 
WRITE (4,40) 
WRITE (5,50) 



20 


FORMAT 


(IX, F8. 


30 


FORMAT 


(IX, F8. 


40 


FORMAT 


(IX, F8. 


50 


FORMAT 


(IX, F8. 


15 


CONTINUE 



TS ( I ), TB ( I ), HF ( I ) , NUSEX ( I ) , RAYD ( I ) , NUSC ( I ) 

TF ( I ) , BETA ( I ) , DENF ( I ) , DENV ( I ) , K INV ( I ) , CP ( I ) 

K ( I ) , HFG ( I ) , ALF ( I ) , PR ( I ) , SURT ( I ) 

Y ( I ) , X ( I ) , CSF ( I ) 

2, 2X, F9. 1, 2X, F12. 2, 2X, F8. 3, 2X, F12. 7, 2X, F15. 3 ) 
4, 2X, F9. 4, 2X, F9. 1, 2X, F8. 4, 2X, E10. 3, 2X, F15. 3 ) 
4, 2X, F14. 4, 2X, E10. 3, 2X, F14. 4, 2X, F9. 5 ) 
4,2X,F14.4,2X,E10.3 ) 
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